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pi • Abstract. We investigate the action of time dependent detunings upon the excitation 

^ . inversion of a Cooper pair box interacting with a nanomechanical resonator. The 

method employs the Jaynes-Cummings model with damping, assuming different decay 
• — 

Co ' rates of the Cooper pair box and various fixed and t-dependent detunings. It is shown 

that while the presence of damping plus constant detunings destroy the collapse/revival 
^ ^ effects, convenient choices of time dependent detunings allow one to reconstruct such 

^ . events in a perfect way. It is also shown that the mean excitation of the nanomechanical 

O I resonator is more robust against damping of the Cooper pair box for convenient values 

of t-dependent detunings. 



Q> ■ PACS numbers: 03.65 -w, 03.65 Yz, 85.85. +j 

cn 

' 1. Introduction 

\6 

^ ■ A popular and exactly soluble model in quantum optics is the Jaynes-Cumming model 
(JCM). It describes the interaction of a two- level atom with a single-mode of the 
electromagnetic field [H El |3l HI [5l EJ [7] . Over the last two decades various extensions 
^ . of the ordinary JCM have been used in various directions, e.g., as adapted to: (i) the 
study of interaction of a three-level atom with a two- mode squeezed vacuum [8]; (ii) 
the study of atom-field interaction in the presence of a cavity damping [9]; (iii) the 
same as in (i), including an additional (nonlinear) Kerr medium [10]; (iv) the two- 
level atoms inside a cavity acted upon by an external field control [H]; (v) study of the 
nonlinear dynamical evolution of a driven two-photon Jaynes-Cummings model [13]; (vi) 
the study of a generalized Jaynes-Cummings models, including dissipation [TU [15l [16] 
and multiphoton interactions [13 [18]; etc. In all these cases, with interest either on the 
field or on atomic properties, the theoretical approach traditionally assumes the atom- 
field coupling as a constant parameter. Comparatively, the number of works in the 
literature is very small when one considers such coupling and the atomic frequency as 
time dependent parameters [191 1201 [211 1221 [231 [12] , including time dependent amplitudes 
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|24] . However, this scenario is also relevant; for example, the state of two qubits 
(qubits stand for quantum bits) with a desired degree of entanglement can be generated 
via a time dependent atom- field coupling [25]. Actually, such coupling can modify 
the dynamical properties of the atom and the field, with transitions that involve a 
large number of photons [26]. In general, these studies are simplified by neglecting the 
atomic decay from an excited level. Theoretical treatments taking into account this 
complication of the real world may employ a modified JCM. In these expected, 
the state describing the system decoheres, since the presence of dissipation destroys the 
state of a system as time flows. 

In the present work we extend what we have learned from the JCM applied to the 
atom-field interaction to investigate a more advantageous system, from the experimental 
viewpoint (faster response, better controllability, and useful scalability for quantum 
computation [27] ) by considering a nanomechanical resonator (NR) interacting with a 
Cooper pair box (CPB) [T71 [151 12H]- Such nanodevice has been explored extensively in 
the hterature, e.g., to investigate: (i) quantum nondemolition measurements [291 130] . 
(ii) decoherence of nonclassical states, as Fock states and superposition states describing 
mesoscopic systems [3ll [32], etc. The fast advance in the technique of fabrication in 
nanotechnology has implied great interest in the study of the NR system in view of 
its potential applications, as a sensor - to be used in biology, astronomy, quantum 
computation, and in quantum information [33l [311 [35] , to implement the quantum qubit 
[36] and in the production of nonclassical states, as Fock state [371 [38], Schrodinger's cat 
state [39], squeezed states [IQ], clusters states [H], etc. In particular, when accompanied 
by superconducting charge qubits, the NR has been used to prepare entangled states 
[42] . Zhou et al.[40] have proposed a scheme to prepare squeezed states using a NR 
coupled to a CPB qubit; in this proposal the NR-CPB coupling is under an external 
control while the connection between these two subsystems plays an important role in 
quantum computation. Such a control is achieved via convenient change of the system 
parameters, which can set "on" and "off" the interaction between the NR and the CPB, 
on demand. 

In this report we will investigate the CPB excitation inversion, its control, and the 
average photon number in the NR. We will consider dissipation in the CPB due to a 
decay rate from excited to ground states. We will also verify in which way the time 
dependence of the CPB-NR coupling modifies these two properties. To this end we must 
solve the time evolution of the whole CPB-NR system, via the approach presented in 
the following Section. 

2. Model Hamiltonian for the CPB-NR system 

A Josephson charge qubit system has been used to couple with a NR. Here we study a 
modified model where a CPB is coupled to a NR, as shown in Fig. [1] below. The scheme 
is inspired by the works of Jie-Qiao Liao et al. [36] and Zhou et al. [40j where we have 
substituted each Josephson junction by two of them. This creates a new configuration 
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that includes a third loop. A superconducting CPB charge qubit is adjusted via a 
voltage Vi at the system input and a capacitance Cg. We want that the scheme attains 
an efficient tunneling effect for the Josephson energy. In Fig. [T]we observe three loops: 
one great loop between two small ones. This makes it easier controlling the external 
parameters of the system since the control mechanism includes the input voltage Vi plus 
three external fluxes $r and $e(^)- In this way one can induce small neighboring 
loops. The great loop contains a NR which is modeled as a harmonic oscillator with a 
high-Q mode of frequency Vt and its effective area in the center of the apparatus changes 
as the NR oscillates, which creates an external flux $e(^) that provides the CPB-NR 
coupling. 

In pursuing the quantum behavior of a macro scale object the nano scale mechanical 
resonator plays an important role. At sufficiently low temperature the zero-point 
fluctuation of the NR will be comparable to its thermal Brownian motion. The detection 
of zero-point fluctuations of the NR can give a direct test of Heisenberg's uncertainty 
principle. With a sensitivity up to 10 times the amplitude of the zero-point fluctuation, 
LaHaye et al |13] have experimentally detected the vibrations of a 20 MHz mechanical 
beam of tens of micrometres size. For a 20MHz mechanical resonator its temperature 
must be cooled below ImK to suppress the thermal fluctuation. For a GHz mechanical 
resonator a temperature of 50mK is sufficient to effectively freeze out its thermal 
fluctuation and let it enter the quantum regime. This temperature is already attainable 
in dilution refrigerators. 

In this work we will assume the four Josephson junctions being identical, with the 
same Josephson energy the same being assumed for the external fluxes and 
i.e., with same magnitude but opposite sign: $l = — $r = ^x- This interaction actually 
couples the two subsystems. Together with the free Hamiltonian of flux qubit and NR, 
the Hamiltonian of the whole system reads 

H = na^a + I^Ng -]^(^z~ cos cos a„ (1) 

where a^(a) is the creation (annihilation) operator for the NR excitation, with frequency 
Vt and mass m; E^j and E^. are respectively the energy of each Josephson junction and 
the charge energy of a single electron; Cg and Cj are the input capacitance and the 
capacitance of each Josephson tunnel, respectively. $o = ^/2e is the quantum flux and 
Ng = CgVi/2e is the charge number in the input with the input voltage Vi. We have 
used the Pauli matrices to describe our system operators, where the states |0) and |1) 
represent the number of extra Cooper pairs in the superconducting island. We have: 
a, = |0) (0| - |1) a, = |0) (1| - |1) (0| and Ec = eV {Cg + 4Cj) . 
The magnetic flux can be written as the sum of two terms, 

$^ = $1 + Bix, (2) 

where the ffist term $i is the induced flux, corresponding to the equilibrium position of 
the NR and the second term describes the contribution due to the vibration of the NR; 
B represents the magnetic field created in the loop. We have assumed the displacement 
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Figure 1. Model for the CPB-NMR coupling. 



X described as x = Xo{a^ + a), where xq = \JmVt/2 is the amphtude of the oscillation. 

Substituting the Eq.(l2]) in Eq.([l]) and controlling the flux $i we can adjust cos (^^) = 0, 
and making the approximation irBix/^o << 1 the above Hamiltonian results as (in 
rotating wave approximation), 

H = Qa^a + -Wq^z + Ao(o"+a + a^o-_), (3) 

where the constant coupling Aq = —AEj cos (^^^^ ( ""^g" ) ^'^'^ effective energy 
Uo = 8Ec (^Ng — . An important advantage of this coupling mechanism is its easy 
and convenient controllability. 

Next, we will extend the previous approach to a more general scenario by 
substituting n ^ cu{t) = n + f (t) and Aq ^ A(t) = Aq [1 + / (t) /n] [261 El SS]; in 
addition we assume the presence of a constant decay rate 7 in the CPB; Uo is the 
transition frequency of the CPB and Aq stands for the CPB-NR coupling. a± and cr^ 
are the CPB transition and excitation inversion operators, respectively; they act on the 
Hilbert space of atomic states and satisfy the commutation relations [o"+,(3"_] = 
and [(3"^,o--|-] = ±cr±. As well known, the coupling parameter X{t) is proportional 
to Ju{t)/V (t), where the time dependent quantization volume V (t) takes the form 
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V {t) = Vq/ [1 + f (t) /Q] 1^ 1221 US]- Accordingly, we obtain the new {non hermitian) 
Hamiltonian 

1 7 
H = uj{t)d^a + -ojqct^ + A(t)(o-+a + a^a_) - |1) (1| . (4) 

It is worth remembering that non Hermitian Hamiltonians (NHH) are largely used 
in the literature. To give some few examples we mention: Ref. [U], where the authors 
use a NHH and an algorithm to generalize the conventional theory; Ref. ^8], using 
a NHH to get information about entrance and exit channels; Ref. |19], using non 
Hermitian techniques to study canonical transformations in quantum mechanics; Ref. 
|50j . solving quantum master equations in terms of NHH; Ref. |5Tj, using a new approach 
for NHH to study the spectral density of weak H-bonds involving damping; Ref. |52] , 
studying NHH with real eigenvalues; Ref. [53], using a canonical formulation to study 
dissipative mechanics exhibing complex eigenvalues; Ref. |5l], studying NHH in non 
commutative space, and more recently: Ref. [55], studying the optical realization of 
relativistic NHH; Ref. [22] , studying the evolution of entropy of atom- field interaction; 
Ref. [21], using a damping JC- Model to study entanglement between two atoms, each 
of them lying inside different cavities 



3. Solving the CPB-NR system 

The state that describes this time dependent system can be written in the form 

where |0,?t,) {\l,n)) represents the CPB in its excited state |1) (ground state |0)). 
Taking the CPB initially prepared in its excited state |1) and the NR in a coherent 
states I a), and expanding coherent state component in the Fock's basis, i.e., |a) = 
exj9(-|a|V2) E^^=o(«"/v^)l^), we have |a) = E^=o \n). Assuming the NR and CPB 
decoupled at t = and the initial conditions Co,„ (0) = and I]J^Lo l^i," = 1 we 
may write the Eq. ([5]) as |^ (0)) = E^=o |1, ^) • 

The time dependent Schrodinger equation for this system is 

»«^ij|*W>. (6) 

with the Hamiltonian H given in Eq. (j4j). Substituting Eq.dl]) in Eq.(l6]) we get the 
(coupled) equations of motion for the probabilitity amplitudes Ci,„(t) and Co,n+iit): 

= - ma;(t)Ci,„(t) - '-cooCi,n{t) - t\{t)V^Co,n+i{t) - |Ci,„ 

= - Z(n + l)co{t)Co,n+l{t) + '-UoCo,n+l{t) - lX{t)V^C,,n{t) . 

The numerical solutions of the coefficients Ci,„(t), Co,„,+i(t) furnish the quantum 
dynamical properties of the system, including the CPB-NR entanglement. 

As well known, in the presence of decay rate 7 in the CPB the state of the whole 
CPB-NR system becomes mixed. In this case its description requires the use of the 
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density operator pcN^ which describes the entire system. To obtain the reduced density 
matrix describing the CPB (NR) subsystem we must trace over variables of the NR 
(CPB) subsystem. For example, pnr = TrcpB^PcN)'- 

PNR = T.Z'=o [CiAt)ClAt) + CoAt)ClAt)] \n) {nt\ . (9) 
4. Excitation Inversion of the CPB 

The CPB excitation inversion, here denoted as I{t)cPB, is an important observable of 
two level systems. It is defined as the difference of probabilities of finding the system in 
the excited and ground state, as follows 

mCPB = Er=0 [\CUt)\" - |Co.n+l(t)r] . (10) 

The Eq. ( fTOl) allows one to look at the time evolution of the CPB excitation 
inversion. First, we assume the resonant case (/(t) = 0) for different values of the decay 
rate 7, with a = 5 and Q = uq = 2000Ao and assuming the NR initially in a coherent 
state with the average number of excitations (n) = 25 as shown in Fig. |2j With the 
exceptions of amplitudes, the plots (a), (b) and (c) of Fig. |2]show identical for collapse- 
revival effects: the higher the decay rate the lower is the amplitude of oscillation. In 
presence of detuning, where f{t) = A = const, where A ^ ujq, Q, we see that the 
excitation inversion in Fig. [3t^a) occurs within the interval 30 < Xot < 50, whereas in 
Fig. [3]^b) it occurs in the range 60 < Xot < 75 and in Fig. [3]^c) the excitation inversion 
goes to zero rapidly. Considering the case of variable detuning, f(t) = 77 sin(a;/t), we see 
in the plots (a), (b) and (c) of Fig. |3]that the excitation inversion occurs frequently in 
Fig. Hl^a) and desapears when the parameter rj increases, as shown in the plots (b) and 
(c). Now, even considering the worst results obtained in the off-resonant cases, with 
detuning A = rj = 6OA0, as shown in Fig. [3t^c) and Fig. Hl^c) we see that we can recover 
the collapse-revival effects via the increase of the parameter u/, as shown in the plots 
(a), (b) and (c) of Fig. [51 So, the parameter uf plays a fundamental role in the control 
of collapse and revival effect. 

Fig. |6] shows plots of the mean value of the NR excitations in the presence of CPB 
decay rate for various values of amplitude of oscillations (parameter t]). Plots (d) and 
(h) are for the resonant case: in (d) the decay rate is greater than in (h). Plots (a), (b), 
and (c) are for constant decay rates, with (time independent) detuning that increases 
from (c) — )► (b) — > (a). Finally, plots (e), (f) and (g) are for time dependent detuning, 
with the parameter rj increasing from (e) — )► (f) — t- (g). We note that the three 
plots for time dependent detunings (e), (f) and ( g) are better than those for constant 
detunings (a), (b) and (c): despite all plots are concerned with the same decay rate, 
the first group is more robust against decay. For example, comparing the plots (a) and 
(g): although in (a) the fixed detuning is A = 6OA0 and in (g) maximum detuning is 
^max = Tj = 6OA0 we See that in the last case the average value of the NR decays more 
slowly. One observes in the plot (c) of Fig. |5] that the interval Aot, where the presence 
of the time dependent detuning recovers the collapse- revival effect, coincides with that 
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in the plot (g) of Fig. E] where the mean value of excitation is around 5 times greater 
than in the case of constant detuning (plot (a) of same figure). 
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Figure 2. Time evolution of the Excitation Inversion in the CPB with (n) = 25, 
= cjq = 2000 Aq, for f(t) =0 (resonance) and different values of decay rates 7.' (a) 
7 = O.OlAo, (b) 7 = 0.05Ao and (c) 7 = 0.5Ao. 
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Figure 3. Same as in Fig. \^for 7 ~ O.OSAq and different values of detunings {cf. 
f{t) = A = const.): (a) A = IOAq , (b) A ^ 20Xo, (c) A - 6OA0. 
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Figure 4. Same as in Fig. [2] for 7 = O.OSAq and different time-dependent detunings 
(cf. f{t) = ris\n{ujit)): (a) r/ = IOAq and Luf = O.OSAq, (b) 77 = 20Ao a^d w/ = O.OSAq, 
(c) rj = 60Aq anc? w/ = O.OSAq. 



(a) 



(b) 



(c) 



Figure 5. Same as in Fig. [5] for 7 = O.OSAq and different time-dependent detunings 
(cf. f{t) = rism{ujft)): (a) rj = 60Aq and cu/ = 20Ao, (b) rj = 6OA0 and uj/ = 40Ao, (c) 
77 = 60Aq and ljI = SSAq. 



5. Conclusion 

We have considered a Hamiltonian model that describes a CPB-NR interacting system 
to study the CPB excitation inversion, I{t)cPB, and the average excitation number 
of the NR, {n(t))j^j^. We have also considered the off-resonant case, with various 
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Figure 6. Average number of excitations of the NR, versus time, for {n{0)) — 25, 
fl ^ ujQ — 2OOOA0, for the following cases: (a) off-resonance 7 = 0.05Ao, A — 6OA0, 
(b) off-resonance 7 = O.OSAq, A = 20Ao, (c) off-resonance 7 = O.OSAq, A — IOAq, (d) 
on resonance 7 = O.SAq, (e) off-resonance 7 = O.OSAq, rj = IOAq and ujf = O.OSAq, (f) 
off-resonance 7 = O.OSAq, rj — 20Aq and ujf — 0.05Xo,(g) off-resonance 7 = O.OSAq, 
rj = 60Aq and uf — O.OSAq, (h) on resonance 7 — O.OSAq. 

values of the detuning parameter / (/ = 0; / = A; and / = rjsm{u!ft)) and in the 
presence of CPB decay (about 10 times greater than the (neglected) NR decay). These 
properties are characteristics of the entangled state that describes this coupled system 
for various values of the parameters involved. We have assumed the CPB initially in 
its excited state and the NR initially in a coherent state (see preparation in [I7j). So, 
the following three scenarios were treated: (z) both subsystems in resonance (detuning 
/ = 0); (a) off-resonance, with a constant detuning (/ = A 7^ 0), and (m) with 
a time dependent detuning (/(t) = rjsm^uft)). The results were discussed in the 
previous Section: in resume, concerning the CPB excitation inversion, an interesting 
result emerged: although the presence of a constant detuning destroys the collapse and 
revivals of the excitation inversion, these effects are restituted by the action of convenient 
time dependent detunings - even in the presence of damping in the CPB; concerning the 
NR average excitation number, another interesting result appeared: convenient choices 
of the time dependent detuning f{t) makes the NR subsystem more robust against the 
decay affecting the CPB subsystem. For constant values of detuning, our numerical 
results are similar to others in the literature using a master equation (see, e.g., Ref. 

m)- 

Finally we emphasize that the change in magnetic flux $e (cf. Fig. [1]), due to 
the presence of an external force upon the NR, is the responsible for controlling the 
parameters uj(t) and X{t). 
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